Low-energy X-ray emission 
from the Abell 2199 cluster of galaxies 



Jelle S. Kaastra 

SRON National Institute for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands 

j .kaastra@sron.nl 

Richard Lieu 

Department of Physics, University of Alabama, Huntsville, AL 35899 

Johan A.M. Bleeker and Rolf Mewe 
SRON National Institute for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands 

and 

Sergio Colafrancesco 

Osservatorio Astronomico di Roma, Via dell'Osservatorio 2, T00040 Monteporzio, Italy 



ABSTRACT 

In a recent Letter, Berghofer & Bowyer rediscussed the analysis of BeppoSAX LECS data of 
the cluster of galaxies Abell 2199 as presented by Kaastra et al., in particular the detection of 
a soft X-ray excess. Berghofer & Bowyer stated that their analysis method is better suited and 
does not show evidence for a soft X-ray excess. Here we find it necessary to publish a rebuttal, 
because it can be demonstrated that the method used by Berghofer & Bowyer is oversimplified, 
leading to an erroneous result. As a consequence, their statement that our initial analysis is 
incorrect is invalid and the detection of a soft X-ray excess in A 2199 is still confirmed. 

Subject headings: galaxies: clusters: general — galaxies: clusters: individual (Abell 2199) — X-rays: 
galaxies 



1. Introduction 

The X-ray spectrum of the cluster of galaxies 
Abell 2199 has been studied by many instruments. 
The detection of a soft X-ray excess in this cluster 
was first claimed by Bowyer et al. (1998), based 
upon Extreme Ultraviolet Explorer (EUVE) data, 
although they do not represent any analysis work. 
Kaastra et al. (1999) analyzed the BeppoSAX data 
of this cluster of galaxies and found evidence of 
both a soft and a hard X-ray excess at radii larger 
than 300 kpc. This last analysis was based on spa- 
tially resolved spectroscopy with data from the 
BeppoSAX, EUVE and ROSAT missions. 



However, in a recent Letter, Berghofer and 
Bowyer (2002), (hereafter BB) made categorical 
statements that the analysis of Kaastra et al. is 
flawed. BB say explicitly that "Unfortunately, the 
telescope sensitivity profile used is likely to be 
incorrect," that "BeppoSAX LECS does not de- 
tect an EUV excess when the data are analyzed 
correctly," that "using a procedure better suited 
to the analysis of extended sources we show that 
there is no excess in Abell 2199," and that "these 
findings appeared to support the (incorrect) find- 
ing of an excess in this cluster using EUVE data." 
In this Letter, we discuss the recent conclusions 
of BB, and we show that BB in fact did not use 
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a different "telescope sensitivity profile" and that 
their conclusions are not based upon improved cal- 
ibration knowledge. Therefore, in this Letter, we 
demonstrate that by pursuing their own highly 
nonstandard approach to BeppoSAX data, the re- 
sults of BB are susceptible to many problems. 

For a clear discussion of this controversy, we 
focus on the data analysis procedures and omit 
the discussion on the physical implications of the 
existence of a soft excess, since this has already 
been presented in our earlier work (Kaastra et al. 
1999). 

2. Summary of the data analysis by Kaas- 
tra et al. 

The data analysis method used previously by 
us has been described by Kaastra et al. (1999). 
Briefly, the cluster was divided into seven con- 
centric annuli, centered around the bright cD 
galaxy, with outer radii of 3', 6', 9', 12', 15', 
18', and 24'. Background-subtracted spectra for 
these regions were obtained from the BeppoSAX 
Low and Medium Energy Concentrator Spectrom- 
eters (LECS and MECS, respectively), High Pres- 
sure Gas Scintillator Proportional Counter, and 
Phoswich Detection System as well as from the 
ROSAT Position Sensitive Proportional Counter 
(PSPC) and the EUVE Deep Survey (DS). The 
background for the BeppoSAX LECS and MECS 
instruments was obtained from the standard 
blank-sky observations, and it was verified that 
the background level during our observation was 
not enhanced due to a higher level of particle ac- 
tivity. 

For the study of the soft excess, the most impor- 
tant instruments are the LECS, PSPC, and DS, 
but the data of all instruments in the entire 0.1- 
100 kcV range were used and fitted simultaneously. 

The instrumental point-spread function (PSF) 
of the LECS instrument is a strong function of 
energy (Parmar et al. 1997). For example, at 
0.28 keV, the energy resolution is 32 % (FWHM) 
and the angular resolution 9.7' (FWHM). Both 
scale approximately with the incoming photon en- 
ergy E as E~°- 5 . Moreover, the wings of the PSF 
are strongly non-Gaussian. Thus, at low energies, 
a significant fraction of the flux generated in a 
given annulus extends to the neighboring annuli, 
so that the annular spectra are coupled and need 



to be fitted simultaneously. 

Therefore, in our analysis, we fitted the spectra 
of all eight annuli and all instruments (26 spectra 
in total) simultaneously. In the response matrices, 
the position- and energy-dependent vignetting fac- 
tors were taken into account as well as the effects 
of the (energy-dependent) overlapping PSF for the 
different annuli. 

The spectral model that we used consists of 
a thermal plasma in collisional ionization equilib- 
rium for each annulus. For the inner regions, the 
possible effects of resonance scattering have also 
been taken into account for the iron Ka complex. 
In addition, a cooling-flow model with partial ab- 
sorption has been included in the central annulus. 

After applying this model to the data, we found 
evidence at the 99.99 % confidence level of the soft 
excess as described by Kaastra et al. (1999). In 
particular, in the 6'-12' range, the 0.1-0.3 keV ex- 
cess luminosity above the thermal component is 
25%. In this region, the subtracted background 
is smaller than 15% of the cluster signal, while 
the large scale variations of the background in this 
region of the sky are less than 10 % of this back- 
ground. 

3. Summary of the data analysis by BB 
3.1. Background subtraction 

BB focus much of their attention on the back- 
ground subtraction. However, as we noted above, 
this is not really a matter of concern for the ob- 
served cluster signal since (1) the background in 
the relevant radial range is relatively small and 
(2) its level during the A 2199 observation was of 
average level and not enhanced. 

BB attempted to divide the background into a 
flat, time-dependent particle background plus an 
X-ray background. Much earlier, Parmar et al. 
(1999a, b) already showed that the particle contri- 
bution in the 0.1-0.5 keV band is small, less than 
13 % of the total, and that it even decreased by 
about 15 % over two years (our A 2199 observation 
is in the middle of this period). Also, a study of 
the MECS background (see F. Fiore, D. Ricci, & 
P. Giommi 1997) 1 indicates that the particle back- 

iBepSDC Report on LECS and MECS Dark Earth Back- 
ground (http://bepposax.gsfc.nasa.gov/bcpposax/soft- 
ware / cookbook /rep_dark_497. html) 
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ground does not vary in time by more than 30 % of 
the total background. The remaining background 
is the cosmic X-ray background. Since the LECS 
only operates during satellite night time, any con- 
tribution to the background from scattered solar 
X-rays is negligible (Parmar et al. 1999b). 

However, BB scaled the remaining non-particle 
background for A 2199 and A 1795 by a factor of 
4.4 and 3.9, respectively, and justified the magni- 
tude of the scaling factor by time-variable scat- 
tered solar X-rays. This is grossly inconsistent 
with the negligible level of solar X-ray photons 
mentioned above. Moreover, BB did not at all ex- 
plain quantitatively how they arrived at this scal- 
ing factor. Scaling is complicated, because the 
cluster fills the entire field of view, and thus the ob- 
served radial X-ray intensity distribution is a com- 
bination of cluster X-rays and cosmic background 
X-rays. The cluster profile is unknown a priori 
(containing potentially both thermal and excess 
X-ray photons). We must therefore assume that 
their undertaking was entirely ad hoc. 

3.2. Radial profiles 

BB questioned the reliability of the low-energy 
calibration of the LECS. The low-energy response 
used by us is based on ray trace simulations. 
This is what is available to the general observer. 
BB quote Parmar et al. (1997) regarding a dis- 
crepancy by a factor of 1.5 at low energies (0.18 
and 0.28 keV) between ray trace simulations and 
ground measurements. However, Parmar et al. 
suggested that the discrepancy may be caused by 
the scattered X-ray photons during the ground 
calibration measurements, and that there is yet 
no convincing explanation, and that in-flight mea- 
surements are needed to resolve this. Since neither 
did BB offer an improved in-flight calibration or 
ray-trace model, they essentially use the same cali- 
bration data as we did and therefore if our analysis 
would contain flaws as a result of this effect, so the 
same will be true of theirs. 

BB attempt to avoid potential calibration prob- 
lems (without solving them) by switching to a ra- 
dial profile analysis. First they noted that at low 
energies, the FWHM of the instrument is large, in 
particular in the 0.1-0.3 keV band. Such a non- 
uniform resolution is undesirable if radial profiles 
in different energy bands are compared. There- 
fore, they convolved the radial profile in the high 



energy band (0.5-2.2 keV) with a Gaussian of 9.7' 
FWHM. Using a spectral model for the cluster, 
they then compared this scaled, convolved 0.5- 
2.2 keV profile with the observed, unconvolved 
0.1-0.3 keV band. The comparison shows no soft 
excess, even a small soft X-ray deficit in the cen- 
ter, and this leads BB to the conclusion that our 
analysis is wrong. 

There are several reasons why this simplified 
approach by BB is unacceptable. Specifically: 

1. Convolving the 0.5-2.2 keV image with a 
Gaussian does not render it compatible with the 
resolution of the 0.1-0.3 keV image. This is mainly 
due to the strong, non-Gaussian tails of the instru- 
mental PSF, as we show in the next section. 

2. Degradation of the images by smoothing de- 
stroys essential information. 

3. The radial profile in the 0.1-0.3 keV band 
is also significantly affected by the low energy re- 
sponse of higher energy photons (32 % FWHM at 
0.28 keV). Thus, the radial profile in this band 
is sensitive to spectral variations as a function of 
radius. 

4. The effective vignetting corrections to be 
made are dependent on both energy and position, 
and these have been apparently neglected by BB. 

5. The 0.5-2.2 keV band contains the Fe-L com- 
plex which can be quite strong in moderately cool 
clusters such as A 2199. In particular, the strength 
of this complex depends on both the amount of 
central cool gas as well as metallicity. Since both 
components vary strongly with position, this bi- 
ases the 0.5-2.2 keV flux. 

In particular, item 1 is very important for the 
analysis of the soft X-ray excess, as we show below. 

4. Radial profiles 

In order to understand what really happens in 
the analysis of BB, we present here some simu- 
lated radial profiles. First, we generated a clus- 
ter emission profile, for which we have chosen for 
demonstration purposes a simple /3-model, with 
core radius r c of 2.6' and (i = 2/3, parameters 
that approximately describe the structure of the 
A 2199 cluster of galaxies. This image is then con- 
volved with the monochromatic instrumental PSF 
of the LECS. A sufficiently accurate parameteriza- 
tion for the present purpose, based on the publicly 
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available calibration files, is given by 

F(r) = l-[l + (x/af}-\ (1) 

where F(r) is the encircled energy fraction of 
the instrument. For energies of 0.19, 0.277 and 
0.93 keV the scale height a is 16.27, 14.70 and 
5.30', and b is 10.57, 13.32 and 5.12, respec- 
tively. Finally, we then convolve the high-energy 
image with a Gaussian for which F(r) = 1 — 
exp(— r 2 /2a 2 ), with a = 4.12', corresponding to 
the FWHM of 9.7' as used by BB. 

In Figure 1, we show the results of our analysis. 
In all but one case, we have chosen an effective 
energy of 0.19 keV for the low energy band (around 
the center of the soft 0.1-0.3 keV band of BB), and 
0.93 keV for the high energy band (around the 
effective center of the 0.5-2.2 keV band of BB). In 
all cases, we plotted the ratio of the soft to hard 
band, for an equal number of photons. 

Case A in Figure 1 shows the profile ratio 
for an isothermal cluster, without convolving the 
high-energy band with a Gaussian, thus reflecting 
purely the instrumental PSF. The strong drop in 
the center demonstrates the effect of the broader 
PSF at low energies. It is evident that even for a 
cluster without spectral variations, the radial pro- 
file ratio varies strongly and is nowhere equal to 
unity. 

Case B in Figure 1 is similar to case A, but now 
we convolved the high energy band with a Gaus- 
sian as BB did. Again, for an isothermal cluster 
the ratio varies strongly as a function of r, but now 
with opposite signatures compared with the pre- 
vious case. Apparently, the hard band has been 
smoothed too much. Adopting the approach of 
BB, a soft flux deficit around 10' would indeed 
have been inferred, despite the fact that the clus- 
ter is isothermal. Taking one step further, if BB 
would have corrected their profile with a curve like 
case B, they would have found the soft excess re- 
ported by us! 

Case C in Figure 1 is similar to case B (i.e., 
an isothermal cluster with the BB approach of 
smoothing the high-energy band), but now for 
a low-energy band of 0.277 keV instead of 0.19 
keV. The oversmoothing of the high-energy band 
is more evident here. The comparison of curves 
B and C shows that if radial profiles are used for 
this kind of analysis, the spectral/spatial energy 



distribution should be known and modeled appro- 
priately. This has not been done by BB. 

Case D in Figure 1 simulates the effect of the 
presence of a cooling flow or an abundance gra- 
dient. It shows the theoretical hardness ratio in 
case the core radius of the hard component is 0.8 
times smaller than that of the soft component. 
This mimics, e.g., the case when the Fe-L com- 
plex (in the hard band) is centrally concentrated, 
because of either an abundance gradient or a cool- 
ing flow. The profiles have not been convolved 
with the instrument or the Gaussian. It also has 
an ^20 % excess beyond ~6', thereby mimicking 
approximately the soft excess as found by us in 
A 2199. 

Finally, case E in Figure 1 corresponds to the 
case presented in case D, but now the profiles are 
convolved with the instrument and the hard band 
convolved with the Gaussian. A comparison of 
curve E with D shows that the convolved profile 
is completely different from the original profile, 
and comparing curve E with curve B shows that 
the differences due to a different cluster model 
are partially washed out, thus confirming that the 
method of BB tends to destroy information by 
oversmoothing. 

All these effects are completely neglected by 
BB. They simply compare their results with the 
expected - scaled - softness ratio of 1 rather than 
the appropriate curve. That curve can only be ob- 
tained by a full spatial/spectral analysis as, e.g., 
done by our team. But then the procedure of BB is 
obsolete, since our spectral modeling already gave 
the correct answers. 

5. Conclusions 

We have shown in this Letter that the analy- 
sis method of BB for assessing the soft excess in 
A 2199 is inadequate and leads to erroneous con- 
clusions regarding the presence of a soft excess. 

In fact, the method of BB is oversimplified since 
it neglects the intricacies associated with the po- 
sition and energy dependence of the effective area 
as well as the spatial/spectral dependence of the 
LECS PSF. 

As a consequence of their inadequate analysis, 
the results presented by BB are misleading, in the 
sense that a simple method which makes use of ra- 
dial profile ratios is purported to be better than a 
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full, sophisticated spatial/spectral analysis, with- 
out assessing both methods in a controlled exper- 
iment. As we have shown here their method leads 
to unpredictable, erroneous results. 

The discussion of BB is also misleading be- 
cause they claim that the finding of a soft excess 
by Kaastra et al. is due to calibration problems, 
while BB neither established this nor used a better 
LECS calibration. 

Finally we have shown that the method of BB is 
also wrong because of a misunderstanding of the 
BeppoSAX background and corresponding back- 
ground subtraction errors. 

In conclusion, we have shown in this Letter that 
our previous results on the soft X-ray excess in 
A 2199 are still valid within the uncertainties re- 
lated to the spectral and spatial sensitivity of the 
BeppoSAX LECS detector. It is clear that a fur- 
ther investigation of the presence of the soft excess 
in A 2199, confirming or refuting it, can only come 
either from a thorough recalibration of the Bep- 
poSAX instruments or from forthcoming observa- 
tions with other satellites with higher sensitivity, 
such as the XMM-Newton satellite. 
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Fig. 1. — Simulated radial profile ratios as de- 
scribed in the text. 
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